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ABSTRACT 

It is standard in theoretical neutrino astrophysics to use a broken power law ap- 
proximation, based on the Band function, to describe the average photon flux of the 
prompt emission of Gamma-Ray Bursts. We will show that this approximation overes- 
timates the contribution of high energy 7-rays (and underestimates low energy 7-rays). 
As a consequence models that rely on this approximation overestimate neutrino event 
rate by a factor of ~ 2 depending on Earth's column density in the direction of the 
GRB. Furthermore the characteristic energy of neutrinos that trigger a km'^ detector 
is typically 10^^ eV, higher than previously predicted. We also provide a new broken 
power law approximation to the Band function and show that it properly represents the 
photon spectra. 
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1. Introduction 

It is commonly believed that prompt emission by GRBs is due to synchrotron radiation by 
electrons accelerated in internal shocks associated with relativistic jets (with a bulk Lorentz boost 
r of 100-1000). A review of the theore tical and observa tional status of GRBs is beyond the scope 
of this paper. The reader is referred to iMeszarosI (120061). Th e average prompt photon emission is 
often described by fitting a Band function (|Band et allll993|) . High energy neutrino emission by 
GRBs in coincidence with the prompt 7-ray photons has been proposed as a consequence of GRBs 
being a candidate to produce (UHE) ultra high energy (up to ~ 10^° eV) cosmic rays (jWaxman 
19951 : Ivietri|[l995l) . 



Waxman & Bahcalll ( 1997 ) have calculated the diffuse flux due to GRBs. Their work has been 



further improved by calculating neutrino emission for individual bursts (|Guetta et alll2004|) and by 
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perfonning detailed GEA NT4 simulations of proton-photon interactions in the internal shocks 
(IMurase & Nagatakill2006|) . This body of work supports the hypothesis that neutrino detection by 
km^ Cherenkov detectors will probably be in coincidence with a handful of bright GRBs. 

Previous studies have correlated photon emission with neutrino emission by approximating 
the photon spectra with a broken power law. We will show that this approximation overestimates 
the contribution of high energy photons (and conversely underestimates the contribution from low 
energy photons). This has consequences in the number o f events expected b y neutrino telescopes 
such as IceCube (lAchterberg et alll2007cl) and KM3NET (IKappes et alll2007|) and in the character- 
istic energy of the neutrinos detected. 

In section ^we describe the Band function. In section fj3]we follow the standard calculation 
of HE neutrinos emission from GRBs. In section ^we compare the Band function with its broken 
power law approximation. In section ^we calculate neutrino event rates. And in section ^we 
discuss the consequences of the calculation shown here. 



2. Prompt photon spectra 



The time-av eraged prompt GRB photon flux is often described by fitting a Band function 



( Bandetallll993 ): 
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The parameters of the Band function are the amplitude A^, the (asymptotic) low-energy spec- 
tral index a^, the high-energy spectral index P^ and the peak energy Epi^ of the uFi, distribution. 



Typical values for the parameters above are: a^ ^ —1, (3^ 



A^ ^ 0.001-1 photons cm s keV . The break energy e^ and the peak energy Ep]^ are related 
by: 



-2.2, e^ ^ 



30 keV-1000 keV and 



2 + a^ 



E, 
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3. Photo-pion production of neutrinos 



There are many good descriptions of the processes that lead to HE neutrinos from prompt 
GRB ph otons in the l iteratu re. We will not repeat the full derivation here. Instead we refer the 
reader to lGuetta et all (|2004|) and we only present here the major features relevant to this paper. 
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Protons accelerated in internal shocks interact with GRB photons via the process: 

p + 7 ^ A+ -> 7i+[+n] ^12+ + u^^e+ + iy^ + iy^ + iye (3) 

and 

p + 7 -^ A+ ^ 7r° + p . (4) 

Because the photon-proton interaction has to create a A resonance, given a photon energy 
there is a minimum proton energy for which Eqsl3] and |4] can take place. In Earth's reference 

frame: 

1 mi — ml „ 

Correspondingly the neutrinos resulting from Eql3]have a minimum energy: 

1 1 m\ — ml „ 

£„ > J < .,-.. > i^T^^l^^^ (6) 

where < Xp^^+ >~ 1/5 is the fraction of the energy transferred to the charged pion from the 
initial proton energy and the factor of 1/4 arises because on average each one of the four final 
leptons in Eql3]has the same energy. 

It is customary to approximate the Band function fit to the average photon spectra as a broken 
power law (which we label approximation A): 



dN^^ _ { {iMkev) ' ^7 < e. 



^ ' ' Toofc^ j ^7 ^ S 

Supposing that protons have a power law spectrum dNp/dEp ~ E^"^, the neutrino spectrum 
traces Eq|71 

A ({E/elr^ E<el 

-T^ =^A (E/e^^r^ el<E<et , (8) 

[ (E/ey^iE/et)-' E > et 

where a^, = —(3^ — 3, P^ = —a^ — 3 and the neutrino break energy e^ is taken from the minimum 
energy in Eql6l 

" 20(1 + ^)2 4e^ ^^ 

The spectrum is steeper above ef, because of synchrotron energy losses by charged pions. Typically 
el ^ 10^^ eV. Muons also suffer from synchrotron losses, but following usual approximations we 
ignore this effect. 



The neutrino flux normalization Ai, is obtained by supposing that the neutrino fluence (ignor- 
ing synchrotron losses) is proportional to the bolometric photon fluence, S^: 

A, oc / E,-^dE, oc S,= dE, E,--l. (10) 

J E^j^^'^ dEy Jq dEy 

Here the choice of minimum neutrino energy, i^;™™ is unimportant and the maximum neutrino 
energy _£;™^^ can be set so that the maximum proton energy is comparable to the highest energy 
cosmic rays (E"™^"" ~ 10^° eV & E'™^^ ~ 5 • 10^^ eV ). Depending on the spectral indices aj, and 
Pu, the neutrino normalization A,, is independent or a very weak function of E^^^. 

It is possible to calculate the neutrino flux without approximating the Band function. In this 
case the neutrino spectra is: 

Band f {Eu/elr^ E,<el 

aJ e-^-^-!^^W'ilE^-^){^E,le'lf^ t\<E,<el . (11) 



dE, 
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Here we have approximated the flux above the syncrhotron energy to a power law. This approx- 
imation is correct as long as e^ is sufficiently larger than e|^. The normalization, A^, is obtained 
using EqJTOl 



4. Approximations of the Band function 

Figure [H shows the Band function with the parameters: a^ = —1, /5-y = —2, e^ = 300 keV 
and A^ = 0.01 photons cm^^ s^^ keV^^. Also shown is the corresponding approximation A 
from EqjT] The normalization of approximation A has been chosen so that the bolometric fluence 
matches that of the Band function. This is a natural choice for setting the normalization, because 
the photon fluence is used in EqJTOl 

It is quite clear that approximation A overestimates the contribution from high energy photons 
and underestimates the contribution from low energy photons. The reason why this occurs is 
because the break energy of the broken power law has been forced to match that of the Band 
function. 

It may be desirable (e.g. to simplify calculations) to have an alternative broken power law ap- 
proximation. Also we will use this new approximation, B, to illustrate the deficiencies of approx- 
imation A. A better choice for approximating the Band function is to require that the asymptotic 
behavior (for very large and very small E^) of the broken power law matches the Band function. 



while leaving the break energy for the broken power law a free parameter: 

(12) 
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The value of the ejfective break energy e^ is given by the energy at which the two branches of the 
broken power law are equal to each other: 



\100keVJ V 



t7 \ 



lOOkeV 
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2 + a^ lOOkeV 
Which results in an effective break energy that is independent of the spectral indices 



e/37-"7 (13) 



e, = ^, (14) 
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where e is Euler's number. Figure \T\ also shows approximation B with the same parameters as 
before. 



5. Effect of expected number of events 

Given a neutrino spectrum dN^/dEy, the expected number of events in a neutrino telescope 
is: 

.£— . dN 

iVevt = A^ / dEP,{E,- E^'^)S{E,, e)-^, (15) 



M 



where A^ is the muon effective area (1 km^ for IceCube/KM3NET), Pf^^Ey] E'^^^) is the probabil- 
ity of a neutrino of energy E^ to produce a muon with energy equal or greater than the neutrino 
telescope threshold -E™" (we assume 100 GeV) and S{Ey, 6) is Earth's attenuation. The probability 

P^(E,;E--)isgivenby: 

P^{E,- Ef^) = NAaUK) < R,{E.; Ef^) >, (16) 

where < R^{Ey] E^™) > is the average muon range given a neutrino energy E„ and a muon 
threshold -E™™. Earth's attenuation factor is given by: 

5(E„^)=e-"(^)^^"^(^''), (17) 

where z{9) is Earth's column density as a function of angle and a^ is the total //-matter crossection. 
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We have calculated the expecte d number of ev ents for a km^ detector. We ha ve used CTEQS 
for th e neutrino-matter cross-section (ILai et afcOOOt) . We follow the calculation by (|Lipari & Stanev 
199l|) of avera ge muon range. The Earth colu mn density is taken from the Preliminary Earth Ref- 
erence Model (IDziewonski & Andersonlll98l|) . 



Using approximations A and B and the neutrino spectrum derived from the Band function, 
we have studied a GRB with a photon break energy e^ = 300 keV, located at a redshift z=l and 



with Lorentz bulk boost T = 300. We also set the spectral indices to a^ 



-2 and P^ 



-1. We 



have normalized the neutrino fluence of all three spectra to the same (arbitrary) value as described 
in section ^ For this GRB the effective break energy is e^ = 110 keV. The neutrino break energy 
for approximation A is 5.98x10^^ eV and for approximation B it's 1.63x10^^ eV. For all cases 
we have fixed the synchrotron energy break at 10^^ eV. Figure [21 shows the three neutrino spectra. 
Again it is clear that approximation A is inadequate because it overestimates the contribution of 
low energy neutrinos. 

Figure[3]the ratio approximation A to B of expected number of events as a function cos(6') for 
this example GRB. For all GRB locations in the sky we see that approximation A overestimates 
the expected number of events. For steeper angles the ratio is larger, because for approximation B 
the characteristic neutrino energy is higher and therefore Earth's attenuation is higher. 



6. Discussion 



We have shown that the usual choice to describe the photon spectra in the calculation of 
neutrino fluxes from GRBs overestimates the contribution of high energy photons (and therefore 
contribution of low energy neutrinos is overestimated). This results in a higher exepected event rate 
by a factor of ~ 2 for all models that make this assumption. The actual value of the overestimation 
of the expected number of events depends on the matter column depth that neutrinos must cross 
through Earth in the direction of the GRB. Also we have shown that the typical neutrino energy is 
Ri 10^^ eV. The ch aracteristic energy of neutrinos is a factor of e large r than the values obtained b y 
Guetta et all (12004 ) and a factor of 10 than the average value used by IWaxman & Bahcalll (|1997h . 
For back of the envelope calculations we provide a new approximation to the Band function that is 
adequate for the computation of neutrino spectra. 



Kashti & WaxmanI (|2005h have discussed the effects of muon and pion energy losses leading 



to neutrin o flavor flux ratios (at Earth) different than 1:1:1. Previous AMANDA searches for 
neutrinos (JAchterberg et afcOOVall bh. based on Waxman-Bahcall-like assumptions have argued that 
this effect is not important experimentally because the typical energy of the neutrinos detected is 
close to the neutrino break energy. Because, as we have shown, the effective neutrino break energy 



is close to the typical syncrhotron energy break, the effects described by lKashti & Waxmanl (|2005|) 
must indeed be relevant. 
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Fig. 1. — The band function is compared to approximations A and B. The three plots show 
dN^/dE^, EdN.JdE.f and E'^dN^/dE^. The Integral of EdN^/dE^ (fluence) for approxima- 
tions A and B has been normalized to the fluence of the band function. For this example we have 



chosen A^ = 10 ^ photons cm ^ s ^ keV ^, 



a-. 



-l,/3. 



-2 and e^ = 300 keV. For each plot, 



the two vertical lines indicate e^ (right line) and e^ (left line). 



0) 

nj 
u 

V) 

(0 



10-^^ 



TJ 



UJ 

TJ 



CM 

UJ 



10" 











r 


— Band 

— ■ Approx A 


/ 


f^\ 


" 


Approx B 




■ 




/ / i \ 1 


/ / 


\ 


/ / 


\ 


/ / 


\ 


/ / 


\ 


/ / 


\ 


/ / 


\ 


/ / 


\ 


/ / 


\ 


/ / 


\ 


/ / 


\ 


/ / 


\ ■ 


/ / 

■ ■ ■■'■■■! / 1 . . .. 


...1 '■ 1 1 



10* 10^ 10® 10^ 10^ 
E [GeV] 



Fig. 2. — Neutrino spectra corresponding to the Band function, approximation A and B with the 
parameters used in Fig{Tl The normalization A^ was chosen arbitrarily, but for approximations 
A and B the neutrino fluence, ignoring synchrotron radiation, is the same for all three spectra. 
The neutrino break energy is 5.98 x 10^^ eV for the neutrino spectra derived from the Band funci- 
tion and for the one derived from approximation A and 1.62x10^^ eV for the one derived from 
approximation B. The synchrotron energy loss has been arbitrarily set to 10^^ eV. 
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Fig. 3. — Left: Ratio of expected number of events for approximation A to the neutrino spectrum 
derived from the Band function as a function of angle. The parameters chosen for the test GRB 
match those of Fig[T]and|2l For cos(6') = —1 neutrinos go across Earth, while for cos(^) = 
neutrinos glance Earth along the horizon. The structure in the plot is correlated with the density 
profile of Earth as a function of depth, e.g. Earth's core effect is seen between cos(6') -1 and -0.85. 
Right: Ratio of expected number of events for approximation B to the spectrum derived using the 
Band function. 



